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Abstract 

We have investigated the spin-wave spectrum of copper metaborate, CUB2O4, 
by means of inelastic neutron scattering in the commensurate magnetic phase. We 
have found two branches of spin-wave excitations associated with the two magnetic 
sublattices Cu(A) and Cu(B), respectively. In the temperature regime 10K < T < 
21K, where only the Cu(A) magnetic moments are ordered, the interaction between 
the two sublattices is found to be negligible. With this approximation we have 
determined the 'easy plane' exchange parameters of the Cu(A) subsystem within 
standard spin-wave theory. 
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1 Introduction 



Copper metaborate, CuB 2 04, crystallizes with the space group I42d (D^) 
with 12 formula units in the chemical cell. Although the chemical structure 
was already solved 30 years ago by x-ray measurements ([7], see Fig. 1), it was 
only recently that the surprising details of the magnetic structure have been 
revealed with help of various experimental techniques ([l]-[6]). At Tjv=21K 
the magnetic system undergoes a second order phase transition from a para- 
magnetic into a commensurate antiferromagnetic state with a magnetic prop- 
agation vector ko = (0,0,0), i. e. the dimension of the magnetic unit cell 
is identical with the chemical one. Below T*=10K the propagation vector 

— * — * 

changes continuously from ko = (0,0,0) to k = (0,0,0.15) (r. 1. u.) at 
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T=200mK, showing that the magnetic structure becomes incommensurate 
with respect to the chemical lattice along the tetragonal c-axis. Near T* we 
observed higher order harmonics in the diffraction pattern of neutron diffrac- 
tion measurements [4,5], which indicates the existence of a magnetic soliton 
lattice in this compound without the application of an external field. 

The peculiarities of the magnetic structures are due to the existence of two 
nonequivalent types of Cu 2+ ions, labeled Cu(A) and Cu(B) in the follow- 
ing. The Cu(A) ion is located at site 4b (point symmetry S 4 , 00|) and is at 
the center of a square unit formed by four oxygen ions. Cu(B) is at site 8d 
(point symmetry C2,x||, x=0.0815) and is surrounded by six oxygen ions 
located at the vertices of a distorted octahedron. In the commensurate phase 
(T* = 10K<T<Tat=21K), the Cu(A) moments are confined close to the tetrag- 
onal plane and the value of the moment at T=12K is nearly the same as the 
value of a free Cu 2+ spin of 1/jl b [5]. In contrast, the value of the Cu(B) 
magnetic moment reaches only a value of 0.25 hb at T=12K. Hence, the mag- 
netic state in the commensurate phase seems to be mainly governed by the 
exchange interactions of the 'strong' (in the magnetic sense) Cu(A) subsys- 
tem. Exchange interactions inside the second subsystem, as well as between 
the two spin systems are much weaker and become only important below the 
second transition which occurs at T*, leading to the observed changes in the 
magnetic structure. 

In this work we analysed the spectrum of spin excitations in the commen- 
surate phase in order to determine the exchange parameters in the strong 
subsystem. This entailed much simplification of the complex magnetic sys- 
tem in CuB 2 4 , but our assumptions are justified by the excellent agreement 
between experimental data and theory, as we shall demonstrate below. 



2 Experimental details 

The inelastic neutron scattering experiments were performed at the Institut 
Laue Langevin (ILL), Grenoble, on the thermal triple-axis spectrometer IN22. 
The spectrometer was operated in the constant final energy mode with a 
neutron wave vector A;j=2.662 A -1 . Contamination from higher order neutron 
wavelengths has been suppressed by the use of a pyrolitic-graphite placed in 
the scattered beam. In order to increase the neutron intensity, the collimation 
was removed from the beam and the analyzer was horizontally focussed. In 
this configuration, the energy resolution at the elastic position was about 0.8 
meV. A single crystal of CUB2O4 of about 0.5 cm 3 was mounted inside an 4 He 
cryostat of ILL-type with the [q,q,0] and [0,0, q] crystallographic directions 
in the scattering plane. Typical constant q-scans taken with a counting time 
of ~4 minutes/point are shown in Fig.3. Additional scans were taken on the 
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cold neutron triple-axis spectrometers TASP (SINQ) with a final wave-vector 
kf— 1.64 A -1 and IN14 (ILL) with a final wave-vector kf—1.55 A -1 which 
yielded an improved energy resolution (AE ~ 0.20 meV) and showed that the 
high energy branch of the spin- wave excitations in CuB 2 4 goes to Tiuj = 
when q — > 0. In addition, a low energy branch with a dispersion confined below 
hid ~1.5 meV could be identified. This branch has dispersive modes along the 
c-axis only. 



3 Spin-wave theory 



The analysis of the upper (high energy) branches in the spin-wave spectrum 
of copper metaborate can be done in a first approximation by neglecting inter- 
actions between the different subsystems and assuming exchange interactions 
in the 'strong' subsystem Cu(A) only. In the following, we considered the ex- 
change interactions with four neighbour atoms at the positions (0, ±|,— |) 
and (±|,0,|) relativ to the central Cu(A) ion. Interactions with additional 
neighbours can only exist in long cation chains (see Fig.l,[7]) and do not 
influence the quality of the analysis of the spin-wave spectrum as shown be- 
low. For the same reasons we ignored other possible weak interactions such 
as anisotropic exchanges within the basal plane and dipole-dipole interac- 
tions. Measurements of the magnetization [2] , the magnetic resonance [3] and 
neutron investigations of the magnetic structure [4,5] have shown that Cu(A) 
ions form a two-sublattice antiferromagnetic structure with magnetic moments 
placed almost in the tetragonal plane ("easy plane"). The Hamiltonian with 
S=| is thus of the form 

h = j z £ + j xy J2(s*s* + sfs]) + h(J2 + £ s]). (i) 

ij ij i j 

The parameter h plays the role of an external magnetic field (h = gfisH) 
applied in the basal plane h = h y , or describes the Dzyaloshinskii field where 
z is the tetragonal axis of the crystal and J xy > J z . The Hoist ein-Primakoff 
transformation within the spin wave theory approach [8] for each sublattice 
(xiUiZi) and (xjUjZj) is defined as 

sZj = -s + !>;,(>,, 

sh = /f (&£■ + k:) (2) 
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The relations between the crystal coordinate system and the local axes are 
given by 



x = Xi cos a-j/jSina = — Xj cos a — yj sin a 
y = Xi sin a + yi cos a = Xj sin a — yj cos a 

Z = Zij, 

as shown in Fig. 2. Using the above equations, the Hamiltonian up to the 
second power in the quantization operators becomes 

H = Hq + H\ + H 2 . 

The requirement that the linear part of the Hamiltonian, H 1 , vanishes, leads 
to a condition for the orientation of the equilibrium magnetization of each 
sublattice: 

sine = A (3 ) 

where h± = nSJ xy . n is the number of nearest neighbours, H = —nS 2 J xy N 
the energy of the ground state and N the number of spins in one sublattice. 
For the evaluation of the spin-wave spectra, one has to consider the part in 
the Hamiltonian which is quadratic in the quantization operators, i.e. Hi. 
Substituing Eq. 3 in H 2 yields 

h 2 = ME K*h + E b IW) - B E( fe ^ + - cY.( b i b l + 

i j ij ij 

S S 

B = -(J xy cos2a - J z ), C = -(J xy cos2a + J z ). (4) 
The Fourier transform 

kj = — U E b ik,2k esxp(ikf itj ) 



N k 

where A; is a vector in reciprocal space, leads to the momentum representation 
of the Hamiltonian 

H 2 = E(M^i* + bt k b 2k ) ~ 2B(B k b+ k b 2k + h.c.) - 2C(B k bt k b+_ k + h.c.)), 
k 

B k = cxp(-ik ■ c/4) cos(A; • 6/2) + exp(ik ■ c/4) cos(A; • a/2). (5) 

The factor B k represents the /c-dependence of the exchange interactions in the 
Cu(A) subsystem. A separation into independent branches in the excitation 
spectrum is obtained by expressing the operators of each sublattice with help 
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of collective-excitations operators a k and c k through 

b 2k = -L(a k e ir ' k +c k e^), 

with the requirement for the phases that (p k — r] k = V'fe — Ck — 8 an d exp(2i5) = 
B k /B k , respectively. Using this transformation, we obtain 

H 2 = £((/i± " 2B^Bl)ata k + (h ± + 2S V / ^B*)c+c fe - (6) 
k 

-2C B k B* k (a k a_ k - c k C- k + /i.e.)). 

The diagonalization of the square-law part [8] gives two branches of spin waves 
with Eigenenergies 

E a k c = sj (h ± T 2(B + C)y/Btfi)(h ± =F 2(B - C)^*). (7) 



4 Comparison with the experiment 



For the high symmetry directions with k along [q,q,0], [q,q,q] and [0,0, q] re- 
spectively, and in zero magnetic field (i.e. cos 2a = 1), Eq. 7 simplifies to 



k x , k y and /c 2 are —h, ^-k and — /, respectively, where h,k and 1 are the Miller 
indices. In the isotropic case (J xy = J z ), the two spin wave branches of Eq. 8 
are degenerate. A small anisotropy lifts the degeneracy and the two branches 
split in an acoustic and an optic mode (compare Fig.4). We obtain a good 
agreement between Eq. 4 and the experimental data from a least-squares fit 
for the following values of the exchange parameters: 

nSJ xy = 7.7±0.1meV, J xy = 45K, -p- = 0.95 ± 0.01. (9) 

Jxy 

The obtained spin wave spectra are shown in Fig. 8. From our experimental 
data, we have only clear evidence of one dispersion branch (neglecting the low 
energy branch along c in this analysis). However, numerical calculations based 
on the Random Phase Approximation (RPA) method show that, at the zone 
center, the intensity of the optical branch is weak. Hence, due to limitations 
in intensity and resolution, it appears that exchange values given in Eq. 9 
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represents an upper bound of anisotropy in the Cu(A) magnetic sublattice. 
We note also that with these parameters the theoretical curves for the direc- 
tion /c=(0,0,q) lie somewhat below the experimental points (Fig. 4a). This is 
probably due to the fact that we neglected weak interactions betwee Cu(A) 
and Cu(B), as noted above. The interaction within the Cu(B) subsystem itself 
leads to a second phase transition toward the incommensurate magnetic phase 
at T* ~ 10K. 



As noted above, we observed a second excitation branch with dispersion along 
the (0,0, q) direction only. It is confined to energies much lower than the spin- 
wave branch considered here. This branch is probably related to spin excita- 
tions of the Cu(B) subsystem, which is still disordered above T*. An analysis of 
the dispersion and temperature dependences of these excitations is in progress 
and will be published later. 



5 Conclusion 



We have determined the spectrum of magnetic excitations in copper metab- 
orate in the commensurate magnetic phase at T=12K by means of inelastic 
neutron scattering. Two branches of excitations have been observed with dif- 
ferent characteristics. Whereas the low-energy branch is attributed to fluc- 
tuations within the Cu(B) sublattice, which only starts to order below T*, 
the spin-wave analysis of the dispersion of the upper branch allowed us to 
determine the exchange interactions within the Cu(A) sublattice. The preci- 
sion of the spin model introduced here can be estimated by comparing the 
exchange parameters with values calculated from the molecular-field theory. 
This allows us to estimate the transition temperature from the paramagnetic 
into the antiferromagnetic state. In particular, if this transition is caused by 
the interaction inside the Cu(A) subsystem only, the temperature of the phase 
transition T/v should give values for the exchange parameters which are close 
to the results of our spin- wave analysis [9]. The staggered molecular field on 
every spin gn B ( y H e ^Si + He^S?) is obtained for a single pair interaction JSiS 2 
from 

gii B H eM2 ) = ^{S 2 (i))J. (10) 
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As a result, the relation between the Neel temperature and the exchange 
interaction becomes [9] □ 



ijjrfg+i) 

2 3 v ' 

For T/v = 21if we obtain J = 42K. This value is close to the result of the 
spin-wave analysis, which shows that the theoretical analysis of the spin-wave 
spectrum in copper metaborate in the temperature range T* < T < T/v can be 
described in a first approximation as an antiferromagnetic ordered subsystem 
of Cu(A) ions with an anisotropic exchange of 'easy plane' type. Although the 
spins at the Cu(B) site are still disordered above T=10K, interactions between 
the two subsystems influences the excitation spectrum and leads to a shift of 
the spin-wave spectrum, in particular along the z-axis. 

This work is done under partial support from RFBR (grant 01-02-17270). We 
also thank E. Kusmin, M. Popov and A. Pankratz for fruitful discussions and 
L.P Regnault for his help in setting up the experiment on IN22. 
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Fig. 1. The chemical cell of CUB2O4. Cu(A) and Cu(B) positions are represented 
by black and open symbols, respectively. 
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Fig. 2. Relation between the crystal coordinate system and the local axes of the two 
Cu(A) magnetic sublattices. 
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Fig. 3. Typical constant q-scans taken in CUB2O4 at T=12 K at different positions 
in reciprocal space along the (0,0, q) direction. 
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Fig. 4. Spin-wave spectrum of C11B2O4 at T=12K. Circles - experimental data points 
obtained from inelastic neutron scattering, solid and dash lines - result of the the- 
oretical analysis for (a)-(0,0,q), (b)-(q,q,0) and (c)-(q,q,q) directions, respectively. 
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